Acoustoelectric pumping through a balhstic point contact 
in the presence of magnetic fields 
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The acoustoelectric current, J, induced in a ballistic point contact (PC) by a surface acoustic wave 
is calculated in the presence of a perpendicular magnetic field, B. It is found that the dependence 
of the current on the Fermi energy in the terminals is strongly correlated with that of the PC 
conductance: J is small at the conductance plateaus, and is large at the steps. Like the conductance, 
the acoustoelectric current has the same functional behavior as in the absence of the field, but with 
renormalized energy scales, which depend on the strength of the magnetic field, |B|. 
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I. INTRODUCTION 

Acoustoelectric currents induced in a two dimensional 
electron gas (2DEG) by a surface acoustic wave (SAW) 
have been observed in point contacts (PC's) defii, 
in a GaAs/AlGaAs heterostructure by a split gateEl 
This P|Cobk;ni has been analyzed theoretically in several 
papersBB Q. In particular, it has been realizedQ that the 
acoustoelectric current in a ballistic PC does not simply 
represent the drag of the electrons by the SAW, but con- 
stitutes an example of a pummng phenomenon. Pumping 
in conducting nanostructureaS, first considered ipnBef. |9|, 
is recently a subject-oLeonsiderable theoreticalll3"ll3, as 
well as experimentalEjilll interest. The pumping current 
is excited not by applying a voltage bias to the nanos- 
tructure terminals, but rather by periodically and adia- 
batically varying certain properties of the system, i. e. 
certain parameters of the system Hamiltonian. In the 
case of the acoustoelectric current, the periodic pertur- 
bation is formed by the moving potential profile created 
by the acoustic wave, of frequency uja and wave vector q, 



5U{v,t) — A cos{ujot — qr) 



(1) 



As the frequency of the acoustic wave is small, this per- 
turbation can be considered as adiabatic; moreover, the 
potential SU can be presented as SU (r, t) = Ai (t) cos qr-|- 
A2(t)sinqr with Xi{t) = Acosust and X2{t) — As'mujt 
being the 'pumping parameters'. This form of the 
pumping—Geld is very similar to the one suggested by 
Thoulesfllj, in his pioneering discussion of quantized 
charge transfer. 

This paper is devoted to the calculation of the acousto- 
electric current flowing through a ballistic PC subject to 
a perpendicular magnetic field due to a SAW propagat- 
ing along the PC. We employ a model of non-interacting 



electrons. Our main result is that the acoustoelectric cur- 
rent dependence on the Fermi energy in the leads in the 
presence of the magnetic field B has the same functional 
form as for B = 0, with the energy scales (or alterna- 
tively, the length scales) which characterize the PC, being 
renormalized and replaced by field dependent quantities. 
(This fpatU|re is known for the quantum conductance of 
the POM.) 

The derivation of the acoustoelectric current, in the 
presence of a magnetic field, requires quite a lengthy cal- 
culation. We therefore have organized the paper as fol- 
lows. The next section describes the formulation of the 
problem, and gives the final result for the field-dependent 
acoustoelectric current. We discuss there its dependence 
on the Fermi energy of the terminals, on the magnetic 
field and on the wave vector of the SAW. Section III 
includes the formal derivation of a time-dependent cur- 
rent, induced by a weak time-dependent perturbation, 
in a nanostructure of a general geometry, subject to a 
magnetic field, and (a possible) dc bias. The-iprmalism 
is based on time-dependent scattering statesEll. At the 
end of that section we obtain the pumped acoustoelec- 
tric current, and relate its form to iJae formula for the 
pumping current derived by Brouweio. The expression 
for the acoustoelectric current necessitates the knowledge 
of the full scattering matrix of the PC, in the presence 
of the magnetic field. We carry out this calculation for 
a saddle-point potential in Sec. IV, and derive there the 
result reported in Sec. II. 



II. FORMULATION AND RESULTS 

A quite realistic model for the confining potential 
of a PCriiiJnrtiie plane of a 2DEG is the saddle-point 
potentialQEjy 
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Here m is the electron mass, L is the length of the PC and 
d is its width (fi = 1). For L ^ d this potential corre- 
sponds to a waveguide in the x-direction with parabolic 



walls (at 



< 



L) adjusted to horns (at 



> 



L) with 



opening angle d/L. These horns represent the left and 
right terminals at a; = =pc». The electronic states in this 
confining potential are labeled by the energy E and the 
transverse channel number n. The transmission ampli- 
tude of an electron with energy E propagating in the nth 
channel is 



t{£r. 
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(3) 



where is the (dimcnsionless) energy for the longitudi- 
nal propagation, 



with 



£n^iE- En)/S, En = A(n + 1/2), 



1/mLd, A = l/md^, 



(4) 



(5) 



being the two energy scales of the saddle point potential. 
The energy En is the threshold for the propagation in the 
nth channel: modes having a large positive e„ are prop- 
agating {t is close to 1), while modes with large negative 
En are evanescent (t is exponentially small). 

The motion of an electron through a saddle-point po- 
tential, when a constant magnetic field B is applied per- 
pendicularly to thejix, yj-plane, has been considered by 
Fertig and Halperinllj. Remarkably enough, they find the 
same form, Eqs. ^ and (^), for the transmission coeffi- 
cient, where now the channel index n refers to the Landau 
level to which the edge state belongs, and the two energy 
scales (||) are modified, to be Sb and Ab [see Eq. ( |30| ) 
below] . For example, at strong magnetic fields such that 
ujb = |e|B/mc ^ A, then Ab — ujb while Sb — AS/ujb- 
It therefore follows that, (as in the absence of the field), 
the, conductance (at zero temperature, in units of e'^/h) 
isB g = E„*^(£«), with E [see Eq. (|)] being the the 
Fermi energy in the terminals, Ep. When Ab S> 6b the 
conductance as function of i? is a step-like function, with 
plateaus of width Ab and steps of width Sb- Increasing 
the magnetic field makes the intervals between the steps 
larger and the width of the steps smaller. 

In the following sections we calculate the acoustoelec- 
tric current, in the presence of a magnetic field normal 
to the 2DEG plane, and find that it consists of the sum 
over the contributions from the various channels n, 



J= Jo^F(e„,p). 



Here, Jq is the nominal current, 
Jo = 2ewo-j2-, 



(6) 



(7) 



En is defined as in Eq. (|J) with E = Ep , and p gives the 
SAW wave vector in dimcnsionless units. 
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At zero magnetic field, a — {Ld/2y/'^, while at strong 
fields it tends to eB{L/2dy/'^, where £b = (c/|e|B)i/2 is 
the magnetic length. The function F(s,p), given in Eq. 
( ^ ) below, is the same one which appears in the cal- 
culations of the zero-field acoustoelectric curreniu. This 
function is exponentially small well above and below the 
threshold, where |e| ^ 1. Near the threshold, where 
|e| < 1, this function for p » 1 displays damped oscilla- 



tions as function of p with a period Ap 
to zero at p ^ 0, 



1, and tends 



F(e,p) = 27re-'^^t^(e)c(£) erf 



^ p \ 
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p< 1, 



Here 7^ = 2e lnp-|- argr(l/2 — ie), c(e) ~ 1 and a ~ d/L 
is a screening parameter. It is introduced phenomenolog- 
ically in order to model the screeujng of the SAW poten- 
tial in the wide banks of the PCEl. Since the expression 
(^ for the acoustoelectric current i^he same function as 
in the absence of the magnetic fieldB, but with modified, 
field-dependent, argumopts, we may discuss its proper- 
ties borrowing from ouil3 previous results. 

As in the B = case, the modes whose thresholds En 
are far from the Fermi energy Ep do not contribute to 
the current. This is expected for the evanescent modes; 
for the propagating ones this seems to be "counter- 
intuitive" . It means that in a free channel with almost 
no reflection, the acoustoelectric current is zero. As a 
result, the current is large only when the Fermi energy is 
close to one of the thresholds, \Ep — En\ ^ Sb- Hence, 
as in the B = case, the dependence of the acoustoelec- 
tric current on the Fermi energy closely follows that of 
the conductance. Since Sb < S, increasing the magnetic 
field reduces the domains where the current is strong, 
the squeezing factor at strong fields being wb/A. This 
is in agreement with the expectation that pumping is 
quenched by a magnetic field due to the chirality of the 
electron states. Less expected is that, on the other hand, 
the magnetic field enhances the magnitude of the nomi- 
nal current Jq, Eq. (|^), since Sb < S. This means that 
while away from the threshold the current is reduced, the 
field enhances the current at the threshold, the enhance- 
ment factor at strong fields being (wb/A)^. The field 
changes also the oscillation period in q, which is of the 
order a^^ and hence increases with the field. At strong 
fields a/aB=o — djlB- We also note that the acoustoelec- 
tric current in a saddle-point potential PC is invariant 
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with respect to B — s- — B, in agreement with the up-down 
symmetry discussed in Ref. 

The scahng of the characteristic energies with mag- 
netic field is unique to the saddle-point potential and is 
not a generic feature. However, since the saddle-point 
potential is a good approximation to the actual confin- 
ing potential of a PC in a 2DEG, one expects that at 
least the trends predicted by the scaling will be valid for 
a realistic confining potential. 



III. TIME-DEPENDENT RESPONSE 

The acoustoelectric current, driven by a SAW, is an 
example of a dc current induced by an ac perturbation. 
Quite generally, such a current can be obtaiaed using the 
concept of time- dependent scattering state£^. In this sec- 
tion we use this formulation to derive the expression for 
the current flowing in a multi-terminal ballistic nanos- 
tructure of a complicated geometry, which may be sub- 
ject to a dc bias, in the presence of a magnetic field. The 
result is given in Eq. (|9|) below, and is valid up to sec- 
ond order in the strength of the ac perturbation. We then 
specify to the "pumping" conditions, namely, a dc cur- 
rent induced in an unbiased nanostructure, and shoY-that 
the result is a generalization of Brouwer's formulate for 
the pumped charge, for the case of spatially distributed 
pumping parameters. 

In order to use the time-dependent scattering state 
formalism, one writes the current density operator, 
j(r,i) = -(je/2m)*'f(i.,i)V*(r,t) + h.c, with V = 
V — (ie/c)A(r, t) (where A denotes the vector-potential) 
representing the electron fi<jl4 operator \1/ with the time- 
dependent scattering statesEl, 



M'(r,i) 



dE ^ 



{E)xan{E\r,t). 



(10) 



The state Xan{E\r,t) is a solution of the Schrodinger 
equation with a time-dependent Hamiltonian 



n{r,t) =no(r)+SU{r,t), 



(11) 



excited by an electron coming with energy E from chan- 
nel n of terminal a. The solutions are normalized to unit 
incoming flux. Although the time-dependent scattering 
state is labeled with the energy E of the incoming wave, 
it contains components with energies E' ^ E since the 
Hamiltonian is time-dependent. In (^0|), the operator 
aan{E) annihilates an electron of energy E in the nth 
channel of lead a. The thermal average of these opera- 
tors is defined, in general, by the temperature and the 
chemical potential of the terminal connected to the lead, 
i. e., 

(aL(^)aa'„'(£^')> = ^T^^iE - E')Sa,„,a.'n'fa.{E), (12) 



where fa{E) = (e 



(iS-Mo)/fcBT_ 



T) 



is the Fermi distribu- 



next step is to find the average of the current density 
operator. Then, the total current entering terminal /3 is 
given by 



J pit) 



c?r(j(r,i))|r^oo/5, 



(13) 



where the notation r — > oo/3 means that r approaches 
infinity in lead /3, whose cross section is denoted Cfj. 
Therefore, it is sufficient to obtain the states Xan for 
r oo(3. 

In the case at hand. Tig is the Hamiltonian of the PC 
and 6U{r,t) — J du; ex-p{—iujt)SU{r,uj) is the potential 
created by the SAW. A straightforward calculation yields 
that, to second order in the perturbation SU, the time- 
dependent scattering state is 

Xc.n{E\r,t) = e-^""' [xan{E\r) 

+ J dr'duj'e'^"^''SU{r\uj')(^G{E'\r,r') 

dr"dLu"G{E"\r, r")SU{r", u")e~'^"' 



(14) 



where we have denoted E' = E+uj' and E" = E+lj'+uj". 
Here G is the Creen function of the Hamiltonian Tio, and 
Xan{E\r) is the scattering state of Hq, i. e., in the ab- 
sence of the time-dependent perturbation. The asypip- 
totic behavior of the latter at r — > oo/3 was found inc3, 

Xan (i?|r)|r^oo/3 = SapW 

an 

iE\r) 

w+jE\r)Sp„,^UE), (15) 



xG{E'\r",r'))xan{E\r' 



where S/3rn,aniE) is the scattering matrix element from 
an to Pm for the Hamiltonian Tio, and w~ , w'^ , are 
the incoming and the outgoing waves, respectively. The 
asymptotic behavior of the time-dependent scattering 
state (|lj) is derived from this relation and the asymp- 
totic behavior of the Green functioned! , 



G{E\r,r') 



p ^ -^Y.'^X^^^\r)x0n.{E\r'), (16) 



where the "hat" indicates that the magnetic field is in- 
verted. One finds 



iEt 



Sat3W^^{E\r) 



(i?|r,t)|r 
5](^/3™,a„(i?)u;^„(i?|r) 

du.'e-^-''S^;i^jE',E)w+^{E'\r) 
duj'dLu"e-'^'^'+'^"'>* 



tion with the chemical potential iia in terminal a. The 



xS 



(2) 

/3 m 



^^jE",E',E)w+JE"\r))], (17) 



3 



with 



and 



I [ drxf3^{E'\r)5U{v,Lo')xc.n{E\r), 



i J dvdv'xprn{E"\v')5U{v\uj") 
xG{E'\v\Y)5U{Y,J)Xc.n{E\T). 



(18) 



Upon inserting this expression into Eq. (|T^ for the 
current, using Eq. (]T^), one finds that the current is 
given in terms of current density matrix elements for 
waves w^^. Only diagonal matrix elements contribute, 
and using the normalization of w^^j to a unit flux, one 
obtains 



J pit) 



n 

- ^ Ia{E)\Sli,-a,OL7i{E,t)\^'^, 



(19) 



where a time-dependent scattering matrix, S, has been 
introduced. 



{E,t)=Sp [E] 
5(1) / 

+ I dJdJ'e-^(-'+-"^'S^^l(E", E', E). (20) 



du'e-'^''S^P^^^{E\E) 



The result (|19|) gives the current induced in the system 
by a weak time-dependent perturbation to second-order 
in the perturbation. (The first order contribution to "in- 
ternal potential" in terms of the scattering matrix was 
calculated in Ref. [24| ). It hence may serve to obtain a 
general ac response in the presence of a dc bias. In the 
following discussion, however, we confine ourselves to the 
case of a dc response, in the absence of a dc bias. 

To obtain the dc, pumped current we now assume that 
the nanostructure is unbiased, i. e. all terminals have the 
same chemical potential, and select from the general ex- 
pression only the time-independent terms. The resulting 
expression is 



jdc 



dE 
2^ 



dhJ 



(^f{E)~f{E + cj) 



(21) 



Note that the second order contribution to the time de- 
pendent scattering matrix 5*^^-' does not enter this expres- 
sion. To derive it, we have used the following ptpperties 
of the scattering states and the Green functionEj, 



(3m,anXan 



.,l3mXan — XfSm^ 



(22) 



G{E\r2,n)~G*{E\rur2) 
~i'^XaniE\r2)xlniE\ri) 



(23) 



The form (21) of the pumping current can be pre- 
sented as a generalization of the adiabatic pumping for- 
mula given in Ref. |l^. One notes that S'^^\ Eq. (|l|), 
in the adiabatic approximation is just the correction SS 
to the scattering matrix due to the small perturbation 
6U (see also Ref. ^). In terms of scattering states this 
correction is 



(E) 

and consequently. 



Xan 



(24) 



jp — e 



— / dwuj 
27r 



d.f{E)\ 
dE 



j dridr2SU{ri, —uj)6h({r2,uj) 



E- 

amn 



5S. 



f3m,an 



(E) 5Sl^,cn{E) 



5U{v2,uj) 5U{vi,—uj) 



(25) 



It follows that the above derivation gives the pumping 
current for a multi-terminal device, in the presence of a 
magnetic field, which is pumped by a weak modulation 
of the "distributed parameters" 5U (r) . 



IV. THE ACOUSTOELECTRIC CURRENT 

When the time-dependent potential is created by a 
SAW moving along the a:-direction, one has 

5U{r, uj) = Ae'''^5{Lo - wo) + A* e-'''^' 5{lo + ojq). (26) 

Then introducing this into S'^^^ given by Eq. ( p^ one 
finds that at zero temperature, the current pumped by 
the SAW is given by 



amn 

with 



q |2 

/3m, an I 



(27) 



drx/3m(£;|r)e^9^Xa„(i?|r), (28) 



where now E — Ep. 

The evaluation of the integrals ( |2^ ) requires the so- 
lution of the scattering problem of a saddle-point po- 
tential in the presence of a magnetic field. This may 
be accomplished by choosing the symmetric gauge for 
the vector potential, A = (B/2){— y, x, 0}, and chang- 
ing the variables x^px into ^, s, as has been done in 
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Ref. ^ As discussed there, ^ corresponds to the x- 
coordinate of the guiding center, while s describes the 
cyclotron motion about that center. The advantage of 
the {^, s}-representation is that the Hamiltonian, which 
includes the magnetic field, takes again the same form of 
a saddle-point potential. 




^2 1 



but with magnetic field dependent energy scales: 



(5b = 2./72-A2_, Ai 



(29) 



(30) 



That is, in the {^, s}-representation, the variables are 
separated. The parameters of this transformation are 
given by 



X 

Py 



2£ 



e sm ( 



V2£ 



e cos ( 



(31) 



Here, 



£=(2toA)-i/^ tan20 = -tJB/47, 
tanh26'i = -A-/7, tanh26'2 = 7/A+, 



(32) 



where 



1 

— L 

4 
1 

2' 



A±--A±\/72 



1 ^2 



(33) 



The scattering solutions of Hamiltonian (p9|), of energy 
E, are labeled by an, where a = 1,2 correspond to the 
left and right terminals, and and n is the channel in- 
dex. xin{E\Cs) describes an electron entering the PC 
from X — —CO, and X2n{E\^, s) corresponds to an elec- 
tron comipg from x = +00. In the presence of the mag- 
netic fieldt3 xin describes an electron entering from the 
upper left corner (assuming B > 0) of the {a;, y}-plane, 
while X2n belongs to an electron coming from the right 
lower corner. We write the scattering states in the form 



Xan{E\^, s) = C$„(s)xa«(e«|C), a = 1, 2. 



(34) 



Here C is the normalization constant, chosen to be 
real (to be determined below), and <&„ are the eigen- 
functions of the s-part of the Hamiltonian (|30|), i.e. 
harmonic oscillator wave functions, corresponding to 
eigenenergies En = AB(n-|-l/2) and normalized to unity, 
/^c?s$2(s) = 1. Then Xan{£n\S,) are the scattering 
states of the one-dimensional problemil, defined by the 
^-part of the Hamiltonian (|3C|) , with (dimensionless) en- 



ergies Er, 



Xi,2(en|^) = -it(e„)E(-£„,±0, 



(35) 



and E is the complex Weber (parabolic cylinder) func- 
tion, as defined in Ref. |2^. 

The integrals (|2^) are not matrix elements because the 
functions x and x belong to different Hamiltonians, one 
with B and the other with — B. As can be appreciated 
from the above analysis of the saddle-point potential in 
the presence of a magnetic field, the transformation ( |3l| ) 
turns the scattering problem formally into the one in 
the absence of the field, but in the mixed representa- 
tion {momentum, coordinate}. It follows that in order 
to perform the integrals (|2^), it is convenient to convert 
them first into matrix elements, which are then straight- 
forwardly transformed in accordance with ( pT| ) . This can 
be done by using the relations (p2[) to obtain 



E Vn,/3'™'(i?)M|,„,^_(i?), (36) 

/3'm' 



with 



M'p^,aniE) - J drx}UE\r)e^'^''Xc.n{E\r) 

= (/3m|e*«"|a">- (37) 

Hence, the calculation of the acoustoelectric current is 
reduced to (i) obtaining the scattering matrix elements 
(including their phases, which are not derived in Ref. ^|); 
(ii) normalizing the scattering states ( ^ ) to a unit flux; 
and (iii) evaluating the matrix elements (^7|). 

(i) Because of the variable separation in the Hamilto- 
nian (^o|), the scattering matrix docs not mix the modes 
n, 



San.pmiE) — <^rmi>5'a/3 (Cn) , 



(38) 



where Sa/^isn) is the 2x2 scattering matrix for the one- 
dimensional problem defined by the ^-part. The ab- 
solute values of these matrix elements have been de- 
rived in Ref. |l|: |S'i2(e„)| = |5'2i(e„)| = i(e„), and 
|^ii(e„)| = |^22(£„)| = r(£„) = e-^"t(e„). To find tb| 
phases of the matrix elements Sap we use the relatione^ 

E(-£, -0 = -ze-^^E(-e, + tt{e)-'E*i-e, 0, (39) 

which enables us to write X2 in terms of xi and Xi- ^ 
second expression connecting xi with X2 is provided by 
Eq. (^). From these relations, and the fact that xi and 
Xi are linearly independent, we obtain 



S'll(en) = S22{£n) = -?"(£«), 
Sl2{£n) = S2l{£n) = -it{en)- 



(40) 



(ii) In order to deduce the normalization constant 
C in (^4|) we need to obtain the flux transmitted at 
X — +00 due to the scattering state Xin{E\£,., s) coming 
from X = — oo. Utilizing the (brajket) notations, such 
that {xy\lEn) = xin{E\x, y), we write the a;— component 
of the flux at the PC cross-section at 
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Jx{xo) = ^( / i dy{lEn\xy)d{x ~ XQ)vx{xy\lEn) 



n\^{lEn\6{x ~ xn)u^\lEn)j 
Ch\enW(^Jd^ J ds$„(s)E*(-£„,C) 

X [6{X ~ Xo)-ya:]«'n(s)E(-£„, £,) 



Mb J-oo 
xE*(-e,e^)E(-e,ea), a,/3=l,2, 



(48) 



(41) 



This calculation requires the velocity operator, Vx = 
Px/m+{eB/2mc)y in the s} representation. Following 
again Ref. n9| we have 



where we have used Eq. (pa) and the notation ^1^2 = 

Inserting the results (|38|) and (47) into the expression 
for the pumped current induced by the SAW, Eq. ([2^), 
we find 



(49) 



Px 

y 



-e ' cos(p 
-2ie'^^ sine 



1 

V2ie- 



' sm ( 
cos c 



an I3'l3" 

Mp>a{en-,p)M*p„^{£n,p) - {p ^ -p)) ■ 



-id/d£, 
—id/ds 



, (42) 



so that Vx is a linear combination of d/d£, and d/ds. 
Similarly, [see Eq. ([bi])], the operator a; in S{x — xq) is 
a linear combination of f and s. However, at xq —^ +00, 
the term linear in s can be neglected. As a result, the 
terms involving d/ds in Eq. (^) disappear, leading to 



Jx{xo) = CH'{en)2SBa / d^.^^ - xq) 



x3fi 



E*(-e„,e)(-*^)E(-e„,0 



(43) 



in which a, Eq. (^, is the coefficient relating a; to ^ in 
(pT|). Using the asymptotic expression£j of E(a,a;) at 
X +00, enables one to perform the ^-integration, and 
with the requirement Jx{xo) = |>S'i2p, we find 



C 



(2(5b) 



1/2 



(44) 



(iii) The matrix elements ( ^7| ) can be written in the 
s} representation 



/ + CX5 
-00 

X exp[ig(a^ + 6s)]. (45) 
Here b is the coefficient relating x to s in (|3T|), 
1 



A2 - A2\l/2 



mAs A| + SI 



(46) 



It is zero at B=0, increases with B, being of the order of 
the magnetic length £b when ujb — A, and saturates at 
£b when ojb ^ A. Since s ~ 1 and the SAW wavelength 
27r/(7 3> ^Bj the second term in the exponent can be ne- 
glected. Then the s-integration becomes trivial, and the 
resulting expression is 



M'^^^„jE)=S„,,nM0a{en,p) 



(47) 



with 



We parametrise the matrix elements Map using the 
representationcj 

fc±l = + 6-2^^^ T 6"'^% (50) 
where W are the real Weber functions. We then find 



Ml 



■lll = *:^(Kie,p)T^r{s)G{e,p) 

22 ZOb ^ 

M,. = ^(2i7(e,p)±t(e)G(e,p)), 



(51) 



where 



/-I- 00 
d^ (W2(-£, e) + W2(-£, -0) COSpS 
-00 

G(£,p) = Z^" ddw'{-e, -0 - W2(-e,0) sin^, 

J —00 ^ ^ 

H(e,p)- r"rf^W(-£,e)W(-£,-e)cospe (52) 
./ — 00 



Inserting Eqs. (51) and (}40| ) into the current equation, 
(^), we obtain that jf'^ = — •/2'^ = J, where J is given 
by Eq. (^ above, with 

Fis,p)^t^s)Gis,p)H{e,p). (53) 

The integrals G and H where calculated i J. 
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